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D 

0 () 


P 

P' 


rms 


P ij 


r 

Re e 


t 

u.(i = 1, 2, 3) 
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o 

U 

T 
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also used in Cartesian components of fluctuating velocity 
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vector referring to field point position 

Components of field point in Cartesian tensor notation 

vector referring to source point position 

components of source point in Cartesian tensor notation 
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PRESSURE FLUCTUATIONS IN TURBULENT BOUNDARY LAYERS 


SUMMARY 

The equations governing the pressure fluctuations in a turbulent boundary 
layer are derived following the known methods for incompressible flow. The attached 
turbulent boundary layer is discussed, and the hypothesis is advanced that the major 
cause of wall pressure fluctuations is the intermittent eruption of the laminar sub- 
layer observed during flow visualization experiments. This hypothesis agrees with 
all the known features of boundary layer pressure fluctuations and offers a good 
physical framework for the understanding of their actions. 

The separated turbulent boundary layer is also discussed. Compared to the 
attached boundary layer, the equations already derived show how additional sources 
of pressure fluctuations may be present in the separated boundary layer. First, the 
"turbulence-turbulence" interaction terms may be important, and second, it is shown 
how additional pressure fluctuations can arise from the interaction of the turbulence 
with the gradients of velocity in the free stream direction. The exact magnitude of 
these terms will not be known until more detailed experimental information becomes 
available. Finally, it is shown how the conventional turbulence -mean shear inter- 
action gives rise to pressure patterns which are convected at only a fraction of the 
free stream velocity. Since this conclusion contradicts the presently available 
supersonic experimental results, additional mechanisms of pressure fluctuation 
may be occurring in supersonic separated flow. 

SECTION I. INTRODUCTION 


The first work on the estimation of the pressure fluctuations within turbulence 
was due to Heisenberg [Ref. 1] and Batchelor [Ref. 2] . Their work referred to the 
highly idealized model of homogeneous (the same from point to point) , isotropic (having 
no preferred orientation) turbulence. Even for this simple case, considerable manipu- 
lation and approximation were required to derive a result. It was found that the root 
mean square value of the pressure fluctuations in homogeneous isotropic turbulence was 
given by 


rms 


= 0 .58 p Q u 


( 1 ) 


where p Q is the density and ii is the instantaneous velocity fluctuation. Note that the 
magnitude of the velocity fluctuation will be the same for any direction in isotropic tur- 
bulence. 



JLJberoi [Ref. 3] made similar calculations using more detailed experimental 
results. 'These calculations gave an average value of about 
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( 2 ) 


All these calculations have relied on estimating the correlation patterns of the turbulent 
velocity fluctuations, and relating these to the pressure fluctuations via the Navier- 
Stokes equations and all have referred to incompressible flows. Although it is possible, 
in principle, to extend these calculations to include the effects of compressibility, this 
would be a highly complicated task, particularly because of the need to consider re- 
tarded time differences within the flow. 

The results of equations (1) and (2) refer to the pressure fluctuations produced 
by turbulent interactions alone. This source of pressure fluctuations is the "turbulence- 
turbulence" contribution. However, the majority of real flows will contain some mean 
shear, which has a pronounced effect on the pressure fluctuations which occur. This 
second contribution from the "turbulence-mean shear" interaction was first investigated 
by Kraichnan [Ref. 4] . He extended this work [Ref. 5] to cover wall pressure fluctua- 
tions in the turbulent boundary layer. Kraichnan's work has been reviewed and extended 
by Lilley, and Lilley and Hodgson [Ref. 6 and 7] . 

Kraichnan's original calculations in Reference 5 indicated that the ratio of the 
turbulence-turbulence contribution to that of the turbulence-mean shear was 1:32. 
Hodgson has also made calculations [Ref. 8] which give a ratio of 1:20. It would thus 
appear that it is the turbulence-mean shear interaction which gives the major contribu- 
tion to the wall pressure fluctuations, and this result has been generally accepted. How- 
ever, Corcos [Ref. 9] has reported calculations which find the ratio to be only 1:1. 6, 
so that the question cannot yet be considered as completely resolved. 

In the following sections the equations for the boundary layer pressure fluctua- 
tions will first be derived following the known methods, and the various approximations 
required discussed. The contribution of both the turbulence-turbulence and the 
turbulence-mean shear interactions is shown. The attached turbulent boundary layer is 
discussed in detail and a new hypothesis is advanced offering a physical basis for the 
understanding of the pressure fluctuations. The equations are then applied to the sepa- 
rated boundary layer, and it is shown how some of the approximations made for the 
attached boundary layer are no longer valid, so that new sources of pressure fluctuation 
may be expected in the separated turbulent boundary layer. 
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SECTION 1 1, EQUATIONS FOR THE PRESSURE FLUCTUATIONS 


Derivation of the Differential Equation 

In this section the basic equations describing the pressure fluctuations are de- 
rived. This derivation is an amalgam of the methods given in References 4 through 9. 
The analysis begins from the exact equations of aerodynamics which incorporate the 
effects of both viscosity and compressibility, and uses tensor notation with the summa- 
tion convention. 

The equation for the conservation of mass may be written as 


Ip 

at 


+ 


a 

" 5^7 


(p v.) 


= o, 


(3) 


and the equation for conservation of momentum as 


_a 

at 


(p v.) 


•5T (pv i v j > 
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+ XTj (p ij > = 0 - 


(4) 


Both experiment and theory agree that simple non-relaxing fluids show a linear 

dependence of viscous shear on velocity gradient; the aerodynamic stress tensor 

p.. can therefore be written as 
13 

p ij = (p + 'f" pe kl! 5 ij ‘ 2pe ij • 


(5) 


where p is the local static pressure and 

jf±L + JZl l 

e ij 2 { axj ax. ) 
is the "strain" tensor. Als.o, 

a v 

e kk 3x k * 

The derivation of equation (5) is given in detail by Jeffries [Ref. 10] . Differentiating 
equation [ 5] and rearranging gives 
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(6) 


Differentiating equation (4) with respect to Xj /equivalent to taking the divergence) and 
using equation (6) gives 

0. (7) 

Using equation (3) on the first term of equation (7) and interchanging the "dummy" 
suffixes in the last term gives 
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aT 


9(pv ; ) 

3x. 


a 2 a 2 

■ <fv i v J ) + 4r ■ m 


9x .9x . 
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1 9 
3 9x 
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9x.9x. 
1 J 


(P v j v j) 



0. 


( 8 ) 


Equation (8) is still exact within the limitations of equations (5). However, to 
proceed further, it will be necessary to approximate, and from here on the fluid will be 
assumed to be incompressible. With the assumption of incompressibility, equation 
(3) becomes 


9 v. 



9 x. 


i 


0. 


(9) 


When this result is applied to equation (8) , the last term vanishes. It can be seen, 
(therefore, that the viscosity has no direct effect on the pressure, and that even in a 
(compressible flow the effect of viscosity will only enter via the effect of compressibility. 
Thus, the last term in equation (8) may be ignored in all practical cases. This con- 
clusion for local pressure fluctuations parallels that of Lighthill for sound radiation. 
Lighthill [Ref. 11] showed that the inertia terms gave the leading contribution to the 
sound radiation, and in the present case, the inertia terms are again the most important. 
For an incompressible flow, the first term in equation (8) is also zero, although this 
term will certainly require consideration in any treatment of a compressible flow. How- 
ever, this first term will be ignored henceforth. Thus, for incompressible flows equa- 
tion (8) becomes 



9x.9x. 
• J 


(p V.v.) 

' r o I J 



0 . 


(10) 


This is the equation from which the pressure fluctuations may be calculated. 
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Now, 

3^(p v.v.) 
r o i j 


dx.dx. 

' J 
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3x. 
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3v. 

J 

1 3x. 
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+ v. 
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3v. 


3x. 
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( 11 ) 


Using (9), the first term in equation (11) may be seen to be zero, leaving 


3 (p v.v.) 

° 1 J_ 

3x ,3x . 

J 




( 12 ) 


Using equation (9) again, the second term in equation (12) is zero, so that equations 
(10) and (12) show 



3v. 

3v. 

i 

P o 

3x. 

3x. 


(13) 


If the variables are now put equal to the sum of their mean and fluctuating parts 

v. = U. + u. \ 

i i i| 

v. = U. + u. } (14) 

J J J l 

p = P + p' » 

then the relationships for the mean and fluctuating parts may be established. Putting 
relations (14) in equation (9) and taking means gives 


3U. 

^ = o. 


(15) 


since the mean of the fluctuating quantities is zero. Subtracting equation (15) from 
equation (9) shows also 


3u. 

U 1 = °. ( 16 ) 


Now putting relations (14) inequation (13) yields 


5 



2 2 
V P + V P' 



and taking means yields 
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( 17 ) 


(18) 


Subtract in g equation (18) from equation (17) , and noting the identity of the two middle 
terms of the right side of equation (17) , which differ only in their "dummy” subscripts, 
gives 
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Note also that, by virtue of equation (16) , 


a 2 

9 u.u. 
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9x.9xj 
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so that an alternative form of equation (19) is 
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9U. 
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9x. 
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' J 


(19) 


(20) 


( 21 ) 


The right-hand sides of equations (19) and (21) have been split into two parts. The 
first term depends on both the mean shear 91^/ 9xj and the turbulence intensity Uj, 
whereas the second part of the expression depends only on turbulence intensities. These 
two parts correspond to the turbulence-mean shear and turbulence-turbulence contribu- 
tions, respectively. 


Application to the Turbulent Boundary Layer 


Suppose now that equation (21) is applied to the case of the two-dimensional 
turbulent boundary. Ignore, for the time being, the nine component contribution of the 
turbulence-turbulence term. The contribution of the turbulence-mean shear term may 
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I 


be written, in ordinary Cartesian coordinates with x in the free stream direction and y 
normal to the wall, as 




3U 3u 3V 3v 9U 3v 3V 3u 

3x $x 3y 3y 3y 3x 3x 3y 


( 22 ) 


where U, V are the components of mean velocity and u, v are the fluctuating components 
of velocity in the x and y directions, respectively. 

The orders of magnitude, 0() , of these terms can be calculated following the 
classic boundary layer approximation technique. 

Equation (15) shows 



where 6 is the thickness of the boundary layer. Using equation (23) in equation (22) 
gives the terms as 

where the orders of magnitude of each of the fluctuating terms, O(f) , have been taken 
as equal. In the attached boundary layer, the data of Klebanoff [Ref. 12] , reproduced 
in Figure 2, indicate that the u fluctuation is about 50 percent larger than the v fluc- 
tuation, with the magnitude of the lateral fluctuations falling between that of the u and v 
fluctuations. However, the supposition of the equality of O(f) cannot be confirmed or 
denied for the case of a separated boundary layer since no experimental evidence is 
presently available. For the separated case, the supposition can only be regarded as 
a reasonable « priori assumption. 

When the boundary layer thickness (6) is small, it is clear from equation (24) 
that the term is third dominant. Thus, for the attached boundary layer equation (22) may 
be written 
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( 25 ) 



9U 

dy 


9v 

9x 


This equation has been the foundation of a number of attempts to predict attached 
boundary layer pressure fluctuations [Ref. 5 through 9]. However, for the separated 
case, these arguments are not valid, A discussion of the separated boundary layer 
is given in Section 4, but first of all the formal solution to these equations will be 
written. 


Formal Solution of the Equations 


The equations for the pressure fluctuations, from the original of equation (10) 
through its successive forms of equations (13) , (19) ,(21) , (22) , and (25) ,are all of the 
Poisson type. The solution to this equation is well known. If a typical equation is 
written as 

V 2 p = A(x, t) , (26) 

f w 

then the solution in the absence of any boundary is 

pfc, t) = 4- f 4%r d a < 27 ) 

~ 4Tr y x. 1 ~ 

Here x is a vector describing the field point and^ is the "dummy variable" of integration. 
There is an analogy between equation (27) and a solution to the wave equation which 
would be an identical integral, but evaluated at retarded time. This similarity has 
recently been used by Paterson [Ref. 13] in an attempt to produce a simpler mathemati- 
cal and physical model for sound generation by turbulence. However, in the present 
case, we are necessarily interested in the evaluation of the integral at a wall. Here we 
follow Kraichnan [Ref. 5] and include the effects of the wall through a mirror flow model 
of the turbulence field, so that the solution to equation (26) in the presence of a wall 
becomes 


P (x, t) 


4ir 



I A (y, t) | 

rv 



( 28 ) 


where y* gives the value of y reflected in the wall; i.e. , 
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= y 1 

= -y 2 

= y 3 

so that when the field point x is 
becomes 

p( £' f) = -ir /. n 

y 2 >° 

The presence of the wall has caused pressure doubling. Before considering equation 
(29) in more detail, it is useful to return to equation (25) and apply simple similarity 
arguments to give the leading parameters governing the pressure fluctuations. 


ctually on the wall, x - ^ = x - g* , and equation (28) 


4M- d » w . 

Ml 


(29) 


A 

^2 

^3 


SECTION III. THE ATTACHED TURBULENT BOUNDARY LAYER 

Similarity Arguments 


It has been found that the turbulent boundary layer has a velocity profile which 
may be described, at least near the wall (Fig. 1) , by a one parameter family of pro- 
files, the well known logarithmic velocity distribution due to von Karman. This "law 
of the wall" is [ Ref. 14] 



T 


where k and A are experimentally determined constants. U T is the "shear velocity" de- 
fined by 



T 



(31) 


where t q is the wall shear stress. 

Now equation (25) showed the pressure fluctuations to be dependent on the mean 
shear through the flow; this may be found from the velocity profile of equation (30) as 

au u T 

dy k y * (32) 
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In addition, the intensity of the fluctuating velocities through the boundary layer is 
proportional to the momentum transfer through the boundary layer, and thus 

v = a U t , (33) 

where a is a constant of order unity. Using equations (32) and (33) in equation (25) 
shows that 


V 2 p' ~ - p U 2 - (34) 

r ~0 T 

The constant of proportionality in the relationship will depend on the typical eddy size 
and shape, i. e. , on the correlation patterns within the flow. 


Equation (34) establishes an approximate proportionality between pressure fluc- 
tuation and wall shear stress. It has been found experimentally that this proportionality 
gives a satisfactory description of the variations of pressure fluctuation in different 
attached turbulent boundary layers. Lilley has performed detailed calculations [ Ref . 6] 
and finds 


P 


rms 


3. It . 
o 


(35) 


The constant in this equation is found experimentally to be a slowly varying function of 
Reynolds number. The results from a number of experiments are reported by Bull 
[Ref. 15] and generally lie within the range 

2< p /t < 4 . 
rms o 

This simple relationship between pressure fluctuation and local shear is the result of 
the one-parameter family of profiles that can be applied to the equilibrium attached 
turbulent boundary layer. It must be emphasized that any departure from such condi- 
tions will negate the conclusions reached above. For instance, local wall roughness or 
free stream turbulence will affect the results. In particular, there can be no reason to 
apply these results to the very different velocity distributions occurring in a separated 
flow. 

Details of the Pressure Patterns 


The arguments advanced above, although successful, have not given any detailed 
information on the pressure fluctuations. To obtain this information the formal solution 
to equation (25) will be written using equation (29). The pressure fluctuation at the wall 
is thus 

-1 C 2p9U 9 v 

D — / ^ — da . 

P 2ir J n r 9y 9x 
y > 0 
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(36) 



Equation (36) has been written in terms of Cartesian coordinates x, y, z, and r 
is the distance from the field point at the wall. A complete discussion of equation (36) 
would require a consideration of the correlation patterns for the velocity fluctuations. 
However, this refinement will not be considered in this report, and the discussion will 
be limited to a consideration of the local interactions between the mean shear and the 
turbulence. 

To draw some conclusions from equation (36) , the data of Klebanoff have been 
analyzed in detail. These data give the most complete information on velocity fluctua- 
tions presently available, and are thus particularly suitable for the determination of the 
pressure producing mechanisms within the boundary layer. Klebanoff did not make any 
measurements of the pressure fluctuation field, but this is of no particular disadvantage 
in the present case, as a direct comparison of magnitudes is not required. 

Figures 1, 2, 3a, and 3b are from Klebanoff's report [Ref. 12] and show the 
mean boundary layer profile, the velocity fluctuations, and the "dissipation derivatives, " 
respectively. These dissipation derivatives are of particular interest since they also 
represent many of the pressure producing terms of equation (19). The term of particu- 
lar interest is that of 8v/ 9x, as this is the term which interacts with the mean shear in 
equation (36). But this term is one of the smaller derivatives plotted, and it seems 

that the turbulence-turbulence interactions, for instance, of 77^ and ^ , could 

dx ox dzdx 

be important in some cases. The mean shear for this boundary layer has been calculated 
and is shown in Figure 4. This figure has been derived by graphical differentiation of 
Figure 1, and therefore may not be accurate, although it should show the main effects of 
interest. 

Multiplication of 9v/9x from Figures 3a and 3b with 9 U/9y from Figure 4 gives 
a measure of the contribution of each part of the boundary layer to the overall pressure 
fluctuation, and the result is plotted in Figure 5. Clearly, a major source of the pres- 
sure fluctuation arises in the part of the boundary layer closest to the wall, since both 
the mean shear and the velocity fluctuation are increasing rapidly in this region. It 
should be remembered, however, that correlation areas could have a marked effect on 
this conclusion. In Figure 5 the pressure contribution is plotted against local mean 
velocity, as it is thought that this allows a more meaningful interpretation to be made. 
Contributions from the part of the boundary layer below 0. 5 U Q cannot be determined 
since 8v/ 9x is not given there. Measurement of the v component of fluctuating velocity 
requires the use of a cross-wire probe in the hot-wire anemometer, and the physical 
dimensions of the device preclude measurements near the wall. However, in this 
"laminar sublayer" region, the velocity fluctuations (particularly, 9v/9x) may be ex- 
pected to be small. Therefore, the contribution to the pressure fluctuation from well ] 
within this region will also be small. 

Figure 5 is particularly interesting since it shows a major contribution to the 
pressure fluctuation to come from parts of the flow moving with velocities of around 
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0. 6 U Q or less. This contrasts with the generally quoted result that the pressure field 
in the boundary layer is dominated by components moving at 0. 8 U 0 . 

Figure 6, taken from Reference 15, by Bull, offers a solution to this dilemma. 
This figure shows the apparent velocity of convection between two transducers, and 
has been determined from the time delay required for optimum correlation between the 
signals received from each transducer. At large transducer spacings, the apparent 
convection velocity of the overall signal is indeed about 0. 8 U Q , but Figure 6 shows how 
this velocity approaches 0. 53 U 0 at zero spacing. This second figure is in much better 
agreement with the arguments advanced above. 

This variation in apparent convection velocity with transducer spacing has been 
reported by a number of authors [Ref. 16, 17, 18] and is usually explained as being the 
result of the variation in the distance for which various components of the flow are 
coherent. It seems reasonable to suppose that the high frequency small wave length 
components of the flow will lose their coherence in a much shorter distance than the 
low frequency large wave length components. In. addition, it is not unreasonable to sup- 
pose that the small wave length components originate from the slower moving fluid near 
the wall, while the large wave length components are emitted from higher mean velocity 
regions further out from the wall. These two arguments in combination provide an ex- 
planation for the variation in convection velocity with spacing. The relatively rapid loss 
of coherence of the high frequency terms will result in the correlation at large spacings 
being dominated by the faster moving, lower frequency pressure patterns. Additional 
experimental evidence supporting these arguments is presented in References 15, 16, 
and 18. However, some authors have not drawn the inescapable conclusion that local 
pressure fluctuations are dominated by pressure producing phenomena near the wall, 
with a typical convection speed of near 0. 6 U 0 or less. This agrees nicely with the 
arguments from the analysis of the data of Klebanoff put forward earlier. 

The additional experimental evidence supporting these conclusions comes mainly 
from correlation measurements in narrow frequency bands. For this, the transducer 
signals are passed through narrow band filters before correlation. Reasonable con- 
sistency is achieved for high frequency components in the various experimental investi- 
gations, but there is less agreement on the exact effects of the low frequency parts. In 
this discussion the "high" and "low" frequencies may be regarded as being roughly those 
with co <5*/U 0 greater and less than unity, respectively. Bull [Ref. 15] concludes that 
the high frequency components lose coherence after convection for about four wave 
lengths while the low frequency components lose coherence in a way which is not a func- 
tion of wave length. Willmarth and Woolridge [Ref. 16] conclude the loss of coherence 
extends over four to six wave lengths for both high and low frequencies. Both sets of 
results show how the high frequency components are convected at a slow speed, around 
0. 6 U 0 , and the investigations also agree that the convection velocity of the low fre- 
quency components is near 0.8 U Q at large transducer spacings. However, the exact 
effect of the low frequency components at small spacings is not clear. Bull’s results for 



this case seem to lack internal consistency, and Willmarth and Woolridge’s results 
do not a dm it a simple interpretation. This problem is considered further below. 

Possible Flow Mechanisms 

Although the effects of eddy coherence must be of practical significance, the "loss 
of coherence" hypothesis discussed in the last section is not the only mechanism by which 
the pressure fluctuation phenomena may be interpreted. In this section an alternative 
hypothesis is presented which provides a good explanation for many of the phenomena 
observed. Although the arguments are not conclusive, they offer an improved physical 
understanding of the causes of wall pressure fluctuations. 

Suppose that the local pressure fluctuations at any time were dominated by the 
effects of a single eddy structure. Then the variation of convection velocity shown in 
Figure 6 could be simply the result of the outward movement of this eddy as it moved 
downstream. This outward movement would result in its attaining a greater convection 
velocity as it moved downstream. With this idea, Figure 6 may be used directly to 
predict a typical eddy locus from the known boundary layer velocity profile (Fig. 7) . 

The locus represents a weighted statistical average of all the eddies which have 
passed the transducer position. The weighting is provided by the pressure generation 
effects of each eddy, which may be expected to vary with distance from the wall. Thus, 
care should be taken in drawing conclusions from this figure. However, Figure 7 gives 
rise to some interesting speculations about the generation of pressure fluctuations in the 
turbulent boundary layer. It will be recalled that the arguments of the last section 
showed how a major source of pressure fluctuation, particularly of the high frequency 
part of the fluctuation, was located near the edge of the laminar sublayer where the con- 
vection velocities are near 0.6 U Q (Fig. 1). 

In 1956, Einstein and Li [Ref. 19] showed how the laminar sublayer suffered 
from an intermittent disintegration. This conclusion was supported by the work of 
Grant [Ref. 20] and the mechanism has recently been the subject of a thorough investi- 
gation by Runstadler, Kline and Reynolds [Ref. 21] . It is found that the laminar sub- 
layer convolutes itself into eddy structures which erupt away from the wall. Runstadler, 
Kline and Reynolds describe this as "the ejection of momentum deficient fluid from the 
wall" and Grant sees the mechanism as a "stress-relieving motion. " This powerful 
mechanism is at work near to the wall, in just the region where the highest contributions 
to the pressure fluctuations may be expected. It is not unreasonable to suppose, there- 
fore, that a major fraction of the pressure fluctuations at the wall is a direct result of 
this intermittent eruption of the laminar sublayer. 

This hypothesis implies that the pressure-producing eddy is expanding away from 
the wall as it passes downstream, as was concluded above from the convection velocity 
measurements (Fig. 7). Runstadler, Kline and Reynolds [Ref. 21] have made visual 
observations of the eddies resulting from the disintegration of the laminar sublayer and 
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found that the generation is random in space and time. However, they were able to 
define typical paths taken by an eddy after generation as a statistical average of the many 
individual paths observed. This eddy locus is plotted on Figure 7 and may be compared 
with that predicted from the convection velocity measurements. The agreement is 
amazingly good, and provides powerful evidence in support of the hypothesis. The loss 
in agreement at downstream stations may be readily understood as being the result of 
the loss of coherence of the wall eddies. This would result in the convection velocities 
derived from correlation measurements reflecting a larger proportion of the more 
coherent higher speed eddies farther from the wall. But for eddies near to the wall, the 
two plots are almost identical, and one can again infer that the eddies near the wall 
dominate the local pressure fluctuation phenomena. 

However, Figure 7 cannot be regarded as providing conclusive proof of the 
present hypothesis. As has already been mentioned, the locus derived from the fluctua- 
ting pressure convection velocity measurements represents some sort of weighted 
average. In addition, the variation of convection velocity used in the derivation on 
Figure 7 (from Ref. 15) is not identical to that found in other investigations, [Ref. 17 
and 18] . There is also a difference in Reynolds number between the two cases of Figure 
7. Bull's work was accomplished at Reg ~ 20, 000, whereas the work of Runstadler et al, 
had Reg ~ 2, 000. Thus, the agreement between the two cases of Figure 7 could be 
entirely fortuitous. Nevertheless, physical intuition suggests that if such powerful 
"ejection" or "eruption" is taking place near to the wall, then it must be a major factor 
in producing the pressure fluctuations observed. 

The writer feels that this mechanism probably provides all of the high frequency 
part and much of the low frequency part of the pressure fluctuations, and thus implies 
local convection velocities of 0. 6 U Q or less for these fluctuations. The remaining part 
of the low frequency fluctuation probably arises much farther from the wall. It is 
interesting to speculate that this remaining part may occur through an identical eruption 
of the turbulent layer out into the free stream, as part of the known intermittent pro- 
cesses in turbulent flow. The relationship between the eruption of the laminar sublayer 
and the interm ittency observed in the outer parts of the turbulent boundary layer has 
been pointed out in References 20 and 21 and this second phenomenon could well be a 
source of low frequency pressure fluctuations with a convection speed near 0. 8 U Q . It 
is difficult, at this stage, to determine the absolute values of the two postulated sources 
of pressure fluctuation. Examination of the magnitudes of the space time correlation 
curves of Reference 13 suggests that the effect of the inner region through sublayer 
eruption is at least five times stronger than the effects due to the outer region. 

A possible mechanism producing pressure fluctuations, both through sublayer 
eruption and in the outer intermittent region, is eddy acceleration. Eddy accelerations 
have been shown to be a source of noise generation in References 22 and 23. It seems 
probable that these could become important sources of pressure fluctuation in high 
subsonic or supersonic boundary layers. 


SECTION IV. PRESSURE FLUCTUATIONS 
IN A SEPARATED TURBULENT FLOW 


Equation (25) has been the starting point for most of the attempts to predict 
boundary layer pressure fluctuations. However, it cannot be applied to the present 
problem of the separated turbulent boundary layer. First, it is by no means obvious 
that the "turbulence-turbulence" contribution to the pressure (the second part of equation 
19) can be neglected. In a separated flow, high values of turbulence intensity are to be 
expected, and these could contribute substantially to the wall pressure fluctuations. 
However, until definitive measurements have been made in the separated boundary layer, 
it will be difficult to make any final statement on the importance of this contribution. 
Clearly, the turbulence intensities and correlation areas must be expected to be con- 
siderably higher than in the attached boundary layer. There is a general similarity be- 
tween the velocity profile of a separated turbulent boundary layer (Fig. 8) and that of 
a jet. Little work has been done on the pressure fluctuations in a jet, but there has 
been much interest in the parallel problem of the noise radiation. For this case investi- 
gations have shown [Ref. 24] that the turbulence-mean shear interaction is dominant. 

In the absence of other information, it may be assumed, for the present, that the same 
dominance is reflected for the case of pressure fluctuations in separated flow. 

This leaves the turbulence -mean shear interaction of equation (22) as the 
leading source of pressure fluctuations. However, this does not necessarily imply that 
its approximate form given in equation (25) is also valid. Equation (25) has been de- 
rived using the order of magnitude arguments presented in equation (24) . While these 
cause no dissent for the case of an attached turbulent boundary layer, for a separated 
boundary layer, 5 may not be small compared to x. If the separation is taken to occur 
at an angle of 17 degrees, [Ref. 25] then the ratio of 6 to x will be of the order 0. 3. 

The contribution of the first two terms of equation (22) could thus be equal in magnitude 
to that of the third term. This will clearly result in an increase in pressure fluctuation 
and could have a major modifying effect on the correlation and frequency patterns en- 
countered. 

The two regions which may be expected to contribute most of these effects are 
the parts of the flow with maximum positive and negative velocity, where the 9 U/9x 
terms will be highest (Fig. 8). However, little can be said of the exact effects, since 
no information on the structure of the separated turbulent boundary layer is presently 
available. A further complication in the theoretical treatment of these flows is that it 
may not be valid to assume homogeneity in the turbulence field. It is hoped that the pro- 
jected experiments at Wyle Laboratories will help to show the true relative importance 
of the various terms in equation (22) for the separated case. 

It is possible to make some general statements about the effect on the pressure 
fluctuations of the third mean shear term in equation (22) (8 U/3 y) . Figure 8 shows a 
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conjectural velocity distribution for the turbulent separated boundary layer. The mean 
shear becomes large at two separate points in the flow. The first is near the wall. 

Here it may be anticipated that the law of the wall will apply, and that the pressure 
fluctuations from this part of the flow will be proportional to the wall shear stress as 
in equation (34) . This contribution must be expected to be substantially lower than in 
the attached boundary layer since the velocity to which this shear is reacting is the ve- 
locity of the reverse flow, which is only a fraction of that of the free stream. The 
mechanism producing the pressure fluctuations from this region may again be that put 
forward in Section III, i. e. , due to the direct influence of the intermittent eruption of 
the laminar sublayer. 

The second region of high shear is the central part near the "dividing streamline" 
of the separated flow. The contribution of this region may be large. A separated flow 
grows at an angle perhaps ten times that of an attached flow, and this factor of ten is a 
measure of the total momentum transfer that is taking place through the flow. In the 
attached boundary layer, this transfer is a result of only the wall shear stress, but in 
the separated boundary layer, it is the interaction of the forward and reverse flows 
which gives rise to most of the momentum transfer terms. Momentum transfer in this 
second high shear region may thus be expected to be an order of magnitude larger than 
that near the wall in the attached boundary layer case, with proportionate effects on the 
wall pressure fluctuation. 

Lilley and Hodgson [ Ref. 7] have investigated a similar flow with two regions of 
high shear using a wall jet. They also found the major contribution to the pressure 
fluctuations to come from this outer region. A particular feature of this outer region in 
the present separated flow case is that it is moving at a small fraction of the free stream 
velocity. The pressure fluctuation patterns resulting from this region would therefore 
be expected to have a low convection speed. In addition, contributions may be expected 
from the region of high negative velocity through the first two terms of equation (22) , 
and this will again result in the prediction of a low overall convection velocity. 

This conclusion is in opposition to the experimental work of Kistler [ Ref. 26] 
for supersonic separated flows. Kistler's work shows a convection velocity of about 
0. 6 of the local stream velocity. (Note that in his report the lower free stream velocity 
behind the separation shock was not accounted for in the analysis. ) At this stage, it is 
difficult to state where this discrepancy arises. It could be due to the effect of eddies 
in the outer part of the boundary layer entering via the first two terms in equation (22) . 

If this is the case, it will become apparent in the current low speed separated flow ex- 
periments planned at Wyle. On the other hand, the discrepancy could be due to effects 
associated with the supersonic flow, for example, the eddy acceleration effects discussed 
at the end of Section III. It is also possible that Kistler's result is due to moving shocks 
within the separation region. It is clear that the resolution of this discrepancy must form 
a central part of any attempt to explain the pressure fluctuations in supersonic separated 
flows. 
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SECTION V. CONCLUSIONS 


The equations governing the generation of pressure fluctuations in a boundary 
layer have been derived and the various approximations made have been noted. The 
most important approximation is the assumption of incompressibility. The equations 
have been written in a form which allows identification of the contributions of the 
"turbulence - mean shear" and "turbulence - turbulence" to the pressure fluctuation. 

Both of these terms can be significant, but the more important effect is usually the 
turbulence - mean shear interaction. For the pressure fluctuations in a two-dimensional 
boundary layer, the equation takes the form 


V P' = -2 P, 



9 u l 9 V 9 v + 9 U 9 v 
9 x 9 / 9 y 9y9x 


Htr)- («) 


For the attached boundary layer, the most important term will be the third, and 
the data of Klebanoff [Ref. 12] have been analyzed to show the approximate contributions 
of the various parts of the boundary layer. The analysis shows that the major source 
of the pressure fluctuations lies near the edge of the laminar sublayer, where eddy con- 
vection velocities may be expected to be about 0. 6 U Q or less. It is shown how a proper 
interpretation of the results of space-time pressure correlation measurements gives the 
same conclusion. The frequently quoted result of convection velocities near 0. 8 U 0 
refers to the convection speed of the coherent low frequency eddies which dominate the 
correlations at large spacings, but only represent a small part of local pressure fluc- 
tuation. 


It is suggested that the major part of the local wall pressure fluctuations could 
arise from the intermittent eruption of the laminar sublayer observed in flow visualiza- 
tion experiments. This seems physically likely, and the typical eddy path implied by the 
pressure fluctuation measurements following this hypothesis is found to agree with a 
typical eddy path actually observed in experiment, although this agreement may be 
coincidental. A small proportion of the pressure fluctuation, confined to the low fre- 
quency region, is probably due to other parts of the boundary layer, the inner part of 
the intermittent region, for example. A description of the pressure fluctuations as 
being mainly the result of laminar sublayer eruption agrees with all the known-experi- 
mental facts. It provides a good model for a physical understanding of the forces at 
work, and should lead to improved analytical models to describe the pressure fluctua- 
tion phenomena. 

The pressure fluctuations in a separated turbulent boundary layer have been 
discussed in the light of equation (22) . It was shown how all the terms in this equation 
may produce a significant contribution to the pressure fluctuations observed. It was 
also shown that the third term in equation (22) could be expected to yield pressure 
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patterns which were convected at near zero velocity. The other terms would not be 
expected to alter this conclusion, but it is radically different from the experimental re- 
sult of a convection speed near 0. 6 of the local free stream velocity recorded by Kistler 
[ Ref. 26] . This paradox was not resolved, but is thought to represent a fundamental 
problem in the analysis of the pressure fluctuations in separated flow. The possible 
modifying effects of the turbulence - turbulence interaction terms were noted, but again 
no firm conclusions can be drawn until more experimental evidence is available. 
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FIGURE 4. THE MEAN SHEAR GIVEN BY 
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FIGURE 5. CONTRIBUTION TO THE PRESSURE 
FLUCTUATION FROM VARIOUS PARTS OF THE 
BOUNDARY LAYER ESTIMATED FROM 
FIGURES 3a, 3b, AND 4 




FIGURE 6. VARIATION OF CONVECTION VELOCITY WITH SPATIAL 
SEPARATION FROM BULL REFERENCE 15 



FIGURE 7. PREDICTED AND OBSERVED EDDY LOCI AFTER ERUPTION OF 

THE LAMINAR SUBLAYER 
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FIGURE 8. CONJECTURAL VELOCITY DISTRIBUTION FOR THE SEPARATED 
TURBULENT BOUNDARY LAYER SHOWING MAJOR REGIONS 
OF FLUCTUATING PRESSURE GENERATION 


REFERENCES 


1. Heisenberg, W. , "Zur Statistischen Theorie der Turbulenz," Z. Physik, 
Vol. 124 p. 628, 1948. 

2. Batchelor, G. K. , "Pressure Fluctuations in Isotropic Turbulence, " 

Proc. Cam. Phil. Soc. , Vol. 47, p. 359, 1951. 

3. Uberoi, M. S. , "Quadruple Velocity Correlations and Pressure Fluctuations 
in Isotropic Turbulence, ” J. Ae. Sci. , Vol. 20, p. 197, 1953, (also 
NACA TN 3116) . 

4. Kraichnan, R. H. , "Pressure Field Within Homogeneous Anisotropic 
Turbulence," J. A. S. A. , Vol. 28, pp. 64-72, January 1956. 

5. Kraichnan, R. H. , "Pressure Fluctuations in Turbulent Flow Over a Flat 
Plate," J. A. S.A. , Vol. 28, pp. 370-390, May 1956. 

6. Lilley, G. M. , "Pressure Fluctuations in an Incompressible Turbulent 
Boundary Layer," C of A Rep. 133, 1960. 

7. Lilley, G. M. and Hodgson, T. H. , "On Surface Pressure Fluctuations in 
Turbulent Boundary Layers," AGARD Rep. 276, 1960. 

8. Hodgson, T. H. , "Pressure Fluctuations in Shear Flow Turbulence," C of 
A Note 129, 1962. 

9. Corcos, G. M. , "Pressure Fluctuations in Shear Flows," Univ. of Calif. , 
Inst. Res. Berkeley Rep. 183, July 1962. 

10. Jeffries, H. H. , Cartesian Tensors , Cambridge University Press, 1963. 

11. Lighthill,M. J. , "On Sound Generated Aerodynamically I: General Theory, " 
Proc. Roy. Soc. A., Vol 211, pp. 564-587, 1952. 

12. Klebanoff, P. S. , "Characteristics of Turbulence in a Boundary Layer with 
Zero Pressure Gradient," NACA TN3178, NACA Rep. 1247. 

13. Paterson, A. R. , "Incompressible Flows and Sound Radiation, " AGARD 
Rep. 460, April 1963. 

14. Schlichting, H. , Boundary Layer Theory , McGraw-Hill Book Co. , New York, 
4th Ed. 1960. 


24 



REFERENCES (Concluded) 


15. Bull, M. K. , "Properties of the Fluctuating Wall Pressure Field of a 
Turbulent Boundary Layer, " AGARD Rep. 455, April 1963. 

16. Wilimarth, W. W. and Woolridge, C. E. , "Measurements of the Fluctuating 
Pressure at the Wall Beneath a Thick Turbulent Boundary Layer," J. F. M. 

Vol. 14, pp. 187-210, 1962. 

17. Serafini, J. S. , "Wall Pressure Fluctuations and Pressure Velocity Correlations 
in a Turbulent Boundary Layer," NASA TR R-165, December 1963. 

18. Wilimarth, W. W. and Woolridge, C. E. , "Measurements of the Correlation 
Between the Fluctuating Velocities and the Fluctuating Wall Pressure in a 
Thick Turbulent Boundary Layer," AGARD Rep. 456, April 1963. 

19. Einstein, H. A. and Li, H. "The Vicous Sublayer Along a Smooth Boundary," 
Proc. A.S. C.E., Paper 945, 1956. 

20. Grant, H. L. , "The Large Eddies of Turbulent Motion," J. F. M. , Vol. 4, 
p. 149, June 1958. 

21. Runstadler, P. W. , Kline, S. J. , and Reynolds, W. C. , "An Experimental 
Investigation of the Flow Structure of a Turbulent Boundary Layer," AFOSR- 
TN-5421, June 1963. 

22. Lowson, M. V. , "The Sound Field for Singularities in Motion, " Wyle Research 
Report WR 64-10, December 1964. 

23. Ffowcs Williams, J. E. , "The Mechanism of Noise Generation by Supersonic 
Flows," Bolt Beranek and Newman, Rep. 1140. 

24. Lilley, G. M. , "On the Noise From Air Jets," ARC 20,376, 1958. 

25. Lowson, M. V. , "The Aerodynamic Pressure Fluctuations on Space Vehicles," 
Wyle Research Report, WR 64-5, October 1964. 

26. Kistler, A. L. , "Fluctuating Wall Pressure Under a Separated Supersonic Flow," 
J. A.S. A. , Vol. 36, pp. 543-550, March 1964. 


NASA-Langley, 1965 M124 


25 


"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof ” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and technical information considered 
important, complete, and a lasting contribution to existing knowledge. 

TECHNICAL NOTES: Information less broad in scope but nevertheless 

of importance as a contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: Information receiving limited distri- 

bution because of preliminary data, security classification, or other reasons. 

CONTRACTOR REPORTS: Technical information generated in con- 

nection with a NASA contract or grant and released under NASA auspices. 

TECHNICAL TRANSLATIONS: Information published in a foreign 

language considered to merit NASA distribution in English. 

TECHNICAL REPRINTS: Information derived from NASA activities 

and initially published in the form of journal articles. 

SPECIAL PUBLICATIONS: Information derived from or of value to 

NASA activities but not necessarily reporting the results of individual 
NASA-programmed scientific efforts. Publications include conference 
proceedings, monographs, data compilations, handbooks, sourcebooks, 
and special bibliographies. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



